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Abstract 

In his fondamental "Essai sur le probleme des trois corps" (Essay on 
the 3-body problem [Ll| ). Lagrange, well before Jacobi's "reduction of the 
node", carries out the first complete reduction of symetries in this prob- 
lem. Discovering the so-called homographic motions (Euler had treated 
only the colinear case), he shows that these motions necessarily take place 
in a fixed plane, a result which is simple only for the "relative equilibria". 
In order to understand the true nature of this reduction - and of La- 
grange's equations - it is necessary to consider the n-body problem in an 
euclidean space of arbitrary dimension. The actual dimension of the am- 
biant space then appears as a constraint, namely the angular momentum 
bivector's degeneracy. I describe in detail the results obtained in a joint 
paper with Alain Albouy published in french in 1998 ^ AC . : for a non homo- 
thetic homographic motion to exist, it is necessary that the motion takes 
place in an even dimensional space. Two cases are possible: either the 
configuration is "central" (that is a critical point of the potential among 
configurations with a given moment of inertia) and the space where the 
motion takes place is endowed with an hermitian structure, or it is "bal- 
anced" (that is a critical point of the potential among configurations with 
a given inertia spectrum) and the motion is a new type, quasi-periodic, 
of relative equilibrium. Only the first type is of Kepler type and hence 
corresponds to the absolute minimum in Sundman's inequality. When the 
space of motion is odd dimensional, one can look for a substitute to the 
non-existing homographic motions: a candidate is the family of Hip-Hop 
solutions, which are "simple" periodic solutions naturally related to rela- 
tive equilibria through Lyapunov families of quasi-periodic solutions (see 
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|CF) V Finally, some words are said on the bifurcation of periodic central 
relative equilibria to quasi-periodic balanced ones. 

Part of this survey is adapted from JCOf . a course given m Ravello m 
1997. 
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1 From the side of ambiant space to the side of 
the bodies 

This section introduces the disposition space T) and its scalar product fi. This 
leads to distinguishing between the "objective" side of ambiant space, repre- 
sented by an euclidean space {E, e) and the "subjective" side of the bodies, 
represented by the euclidean space (T>,ii). 
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1.1 The shape of N points in an euclidean space 

The equivalent formulae 

^((^i,--- ,^n),{Vi,--- ,Vn)) ■=J2<ri,fj>E^iVj = ^{-^r^j)^ir]j 

'i-j i,j 

define a quadratic form on the hyperplane 

n 

V* :={(^i,...,e„)eM"| $^Ci=0} 

2 = 1 

of M". The vectors fi,--- ,fVi are elements of a finite dimensional euclidean 
space {E, e), whose scalar product and norm are written <, >^ and || ||<;, and the 
positive immbcrs r,j = — rj\\, arc the mutual distances. 

This form defines the configuration of the n points up to a common rigid 
motion (translation and rotation) in E. The hyperplane D* is to be considered 
as the dual of the space of dispositions 

P = M"/(1,--- ,1), 

that is the space of n-tuples of points on a line up to translation. The n{n — l)/2 
numbers rfj are the coordinates of /3 in a natural basis of the space of quadratic 
forms on V*. Considered as a homomorphism from V* to V, /3 takes the form: 

l3 = *xoeoxG Homs{V\V) = Q{V*) = VqV, 

where the isomorphism e : E ^ E* is given by the euclidean structure of E and 
the mapping 

n 

x:V* ^E, a;(a, ■■■ ,^n) = J2^i^i' 

2 = 1 

defines the configuration up to translation (notice that we use the french 'left" 
convention for the transposed mapping *x : E* ^ V** = V). This is nothing but 
the Gram construction. In other words, (3 — x*e. One checks that Ker/3 = Kera; 
and Im /3 = Im*x. A necessary and sufficient condition that the rfj be indeed the 
squares of the mutual distances of n points in some euclidean space E is that 
/3 be positive (Menger [Me] in terms of determinant inequalities, Schoenberg 
[Sch], see also Blumenthal [Bl]). When this condition is not verified, one can 
still interpret the rfj (which one can even replace by not necessarily positive real 
numbers) as the squared mutual distances of n points in some space endowed 
with a non degenerate but not necessarily positive quadratic form. 

1.2 Mcisses as an eucUdean structure on the dispositions 

Given n positive numbers mi,- - ■ , nin, the masses, we define an euclidean struc- 
ture on V by 

n ^ 

\\{xi,---,x„)\\l = '^mi{xi- xcf = rjHmjixi- Xjf, 

i=l l<i<j<n 



3 



where M toi + • • • + m„ is the total mass and xq = (7711X1 + - • • +to„x„)/M is 
the center of mass of the Xi . The associated isomorphism /i : I? — > is defined 

by 

/^(xi, . . . ,a;„) = (mi(a;i - xa), ■ ■ ■ ,?ti„(x„ - xg)), 
'mi m„' 

Notice that, endowing M" with the metric = X]r=i™«^i' one can identify 
the quotient space V with the hyperplane xq = orthogonal to (1, • • • , 1), 
endowed with the induced metric: reducing the translational symmetry amounts 
to fixing the center of mass. All this is in germ in Jacobi; for example, the 
classical Jacobi coordinates amount to choosing a particular orthogonal basis of 
2?*. The quadratic form on E* associated to 

b^xo^o^xe Homs{E*,E) = Q{E*) = EQE 

is called the (dual) inertia form of the configuration. It conveys the same infor- 
mation as the inertia operator! : A'^E* A?E which transforms the instanta- 
neous rotation into the angular momentum (see below) and is defined by 

X{Vl) = bono€-^ +e-^ oflobe [K^E) ®s (A^E). 

The forms /3 and b are closely related : b being on the side of ambient space 
turns with the bodies under an isometry of E, while /3 being on the side of the 
bodies is invariant under isometrics of E. To get a more concrete understanding, 
we transform these forms into the endomorphisms 

B = ^o[3^^o^xoeox: V*^!)* and S = boe = xofj,o*xoe: E^E. 

In an orthonormal basis of E, the endomorphism S is represented by the sym- 
metric inertia matrix 



5': 



I l^k=l^kXik 2^k=l^kXlkX2k 2^k=l'^kXlkXdk\ 

En v-^n 2 

l,^imkX2kXlk l^k=l^kX2k 

\Ek=l^kXdkXlk ■■■ 22k=l^kXik I 



where X\k^ • ■ ■ i x^k £^re the coordinates of Vk—rQ. On the other hand, extended 
to the unique endomorphism of M" sending (toi, • • ■ to 0, B becomes the 

/i~ ^-symmetric intrinsic (or suhjectiv^^ mertia matrix 



B 



I mi||fi - fclle ••• •■• m-i < fi - fG,r*„ - fc >\ 

mi<fi-fG,rj-rG> ■■■ 
\m„ < fn ~ rcfi ~ fa > ■■■ ■■■ Wnl I^Vi - ^^sl I? / 



denomination proposed to the author by the late Adrien Douady 
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The common trace / of the endomorphisms B — ^ o {5 and 5 = 6 o e is the 
moment of inertia of the configuration with respect to its center of mass : 

i—l i<j 

More generally, let volij^...;^ be the volume of the (fc — l)-dimensional paralle- 
lotope generated in E by the vectors r2 — T^iyf^ — fi, ■ ■ ■ , f/j — ri, that is (k — 1)1 
times the volume of the simplex defined by the points ri , ■ ■ ■ ,rk- One checks 
(see [CT]) that 

det{Idv - Mi o /3) = det(Xd£; - A6 o e) = 1 - r;iA + • • • + (-l)""i77„_iA""\ 
where r^u-i = Eii<...<», "^^i ' ' ' "^^A.volL..i, • 

2 The n-body problem and its symmetries 

Given n point masses mi , • • • , to„ located at ri , • • • , rvi in the euclidean space 
(i?, e), the equations of the newtonian n-body problem in a galilean frame are 

They possess the following symmetries: 

1) the galilean symmetry (invariance under a change of galilean frame, which is 
in general dealt with by choosing a galilean frame fixing the center of mass at 
the origin); 

2) the invariance under the action of the orthogonal group of (i?, e); 

3) the scaHng symmetry, which asserts that, given a solution a;(t) = [r\(t)^--- ^rn(t)') 
and any real number A > 0, x\(t) — X^3x{Xt) is also a solution; 

4) the invariance under permutation of equal masses. 

The so-called "reduction" of the n-body problem is directly concerned only 
with the first two but the scaling symmetry is behind the existence of the ho- 
mographic motions with non zero eccentricity. The last one plays an important 
role in the search for "symmetric" action-minimizing (quasi-)periodic solutions 
of the rt-body problem. 

2.1 The reduction of translations: the Wintner-Conley 
matrix 

Since Lagrange |L2| . the equations of the newtonian n-body problem are written 
ruiri = i ^ I,- ■ ■ ,n, = ^ TOimj$(rfj), ^{s) = gs^^. 
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U is the force function, opposite to the potential energy. "Reducing" the trans- 
lations, one can write this equation as the following equality of homomorphisms 
from Vto E* : 



{N) 



eo xo 11 = dU (x) . 



Being invariant under translation, U can be considered as a function x i-^ U{x) 
on Hom(X'* , E) = V^E and its derivative dU {x) as an element of the dual space 
Hom(X'*, £;)*, identified with Hom(D, E*) =T>* ®E* via the bilinear mapping 
{(j), tp) '-^ trace (*^ o tp). The isomorphism x e o x o jj, defines the euclidean 
structure (g) e on }ioui{'D*,E) = 'D^E,m terms of which (iV) becomes 



where V is the gradient. We shall denote by a dot • the scalar product e : 



if X and y are respectively represented by (ri, • • • r„) and (si, • • • s„). For in- 
stance, the moment of inertia I = x ■ x is just the squared norm ||a;||^|g|g of x. 
From the invariance of U under rotations, one can write 



X 



VC/(x), 



n 



X • y = trace {eo x o fio^y) 




U{x) = U{/3). 



From the chain rule one obtains dU = 2eo x o dU, which puts the equations of 
motion in the form 



(iV) 



X = 2x o A. 



where the Wintner-Conley endomorphism 



A = dUoij,-^ : V* 



depends linearly on the masses. Extended to the unique endomorphism of M 
sending (mi, • • • , m„) to 0, A becomes the n x n matrix 



2.2 The reduction of rotations: Lagrange equations 

After reduction of the translations, the phase space of the n-body problem can 
be taken as the tangent space to Hom(I?* , E) and identified with Hom((I?*)^, i?). 
The elements will be written z = {x,y), x for the positions (fi, • • • ,r„), y for 
the velocities {vi, ■ ■ ■ ,Vn)- The space of motion is the image Imz of z and the 
equations of motion are 

(A'') i = y, y = 2xo A. 

By the Gram construction we get (the + sign means symmetric positive) 

^ t f*-xoeox *-xoeoy\ ( B 7 — p\ tt /^^*\2 -t^2n 

£ = *zoeoz^{t t = , . e Hom+ 2?* 2,2?2 . 

yy o eo X y o eo yj \J + P o / 

So, after reduction of both translations and rotations, the state of the system is 
described by four elements /3,7, (5, p of Hom(2?*,I'), whose first three are sym- 
metric (/? and d are moreover positive) and the last antisymmetric. Computing 
£ with the help of (N) we get reduced equations which, under a very compact 
form, coincide with the systems obtained by Lagrange [LI] and Betti [BeJ and 
whose solutions are the relative motions: 

7 = *Ao/3 + /3o^ + (5, 

6 = 2(*A o 7 + 7 o A) - 2(*A o p- po A), 

p = *Ao(3-poA. 

Transforming bilinear forms into endomorphisms of 2?*, we set 

B ~ p,o j3, C = /i0 7, D ^ po 6, R = po p. 

The relative equations become a system of matrix equations 

B = 2C, 

C ^ AB + BA + D, 

D = 2{AC + CA)-2[A, R], 

R = [A,B]. 

2.3 Invariants and first integrals: 1) the energy and the 
Lagrange-Jacobi relation 

The traces /, J, K of the endomorphisms B, C, D can be written 

I = x-x, J ^x-y, K = yy, 
where we recall that the dot means the p ^ e scalar product anV ® E. 
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At the level of traces, what is left of the equations of motion is 

j ^ - ^ K -U, H ^0, 
2 

where H = — U is the total energy, sum of the kinetic energy in a galilean 
frame fixing the center of mass and of the potential energy. The first is the 
Lagrange-Jacobi relation (or viriai relation, see Jacobi [j], Poincare jP] p. 90,91) 
and the second is the conservation of energy. As the Lagrange-Jacobi relation 
can also be written / — AH + 2U, we see that, as I controls the size of the 
system, that is supj^ry, its second derivative I or equivalently the potential 
function U, controls the clustering infi^ rij. 

Remark on the Jacobi-Banachiewitz case. If in the expression of the 
Newton potential we replace $(s) = C/s~2 by = Gs'^ , the Lagrange-Jacobi 
relation becomes / = 2K + AkU = 4iJ + 4(k-|- !)[/, which shows the exceptional 
behaviour of the Jacohi-Banachiewicz potential {k = —1) [Baj IW2| for which 
/ = 2H is constant, which implies the additional first integral G = 2IH — ,P . 
This integral is associated with the scaling symmetry which, in this particular 
case only, is symplectic. 

The Lagrange-Jacobi relation is basic to our understanding of the global be- 
haviour of the solutions of the n-body problem: if we write it J = 2H+2{n+\)U , 
we see that, in the newtonian case or more generally if k > —1, the positivity of 
U implies that J is increasing along any solution whose total energy _ff is > 0. 
The existence of such a Lyapunov function precludes any non trivial recurrence, 
in particular it forbids any periodic motion. It is well known that things are 
much more complicated in negative energy. The basic tool replacing J is, in the 
newtonian case, Sundman's function (see jCO] ) 

s = r^j^ + \c\^)-2r^H, 

where C is the angular momentum bivector which we study in the next section. 

2.4 Invariants and first integrals: 2) the angular momen- 
tum and its associated hermitian structure 

In a galilean frame whose origin is at the center of mass (r^ = "^fc^fc = 

VG = X]fe=i ^kVk = 0), the angular momentum of z = (x, y) G T)^ ® E \s 

n 

C = ^m.ffc A4 e K^E EE HomaiE*,E). 
fe=i 

Considered as an antisymmetric form on E* , that is an antisymmetric homo- 
morphism from E* to E = E** , it can be written 

C = z o oj^ o^z = —X o ^ o *y 4- y o /i o 'x, 
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where : X>2 — t' {V*Y is defined by ujfi{u,v) — (— is the natural 
symplectic structure on T>^ . In other words, 

c = CzYuj^. 

If xik, • • • , Xdk (resp. yife, • • • , j/dfe) are the coordinates of ffc (resp. Vk) in some 
orthonormal basis oi C can be identified with the antisymmetric matrix with 
coeflicients 

N 

Cij = ^ mk{-Xikyjk + XjkUik)- 

k=l 

One readily computes its derivative C — 2x o (— /io*A + Ao/i)o*2: = 0, which 
shows that its invariance follows from the fact that A is /i^^-symmetric. 

The support of the bivector C, that is the image of C G Hom(£'*,i?), is called 
the fixed space. Its dimension, always even, is the rank of C. One finds in |A1| 
the proof by elementary symplectic geometry of the estimates 

rankC < rank£ < -rankC + n — 1, 

which generalize a theorem of Dziobek saying that for three bodies with zero 
angular momentum, the motion necessarily takes place in a fixed plane. 

Hermitian structures. When dimi? = 2 and E is oriented, C can be thought 

of as a real number c. If c ^ 0, the rotation by ±^ according to whether c is 

positive or negative defines on a hermitian structure (i.e. an identification 

of E with the complex plane in which the multiplication by i is an isometry). 

When diuiE = 3, C can be thought of as a vector C orthogonal to the fixed 

plane (if it is different from 0) on which it defines a hermitian structure whose 

multiplication by i is given by the "vectorial product " r -4- A r . 

1^ I 

The same is true in higher dimensions. We first recall the definition of a 
hermitian structure: 

Definition 1 A hermitian structure on the vector space F is a triple 

{e: F ^ F*, J : F F, fl ^ e o J : F ^ F*) 
consisting in compatible euclidean, complex and symplectic structures. 

Definition 2 The hermitian structure (e,J7c,fic) defined on F — ImC by the 

bivector C is defined by 

Jc - [y/-iCoe)^y\coe, ^c^eoJc. 

Of course, in the definition of Jc, we invert \/ —(C o e)^ only on the fixed space 
but the formulas still make sense if we replace F hy E. They define then a 
degenerate (non invertible) hermitian structure, where the action of J7c consists 
in the orthogonal projection on F followed by its action on F. 
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A similar structure (cr, £7, JT") is obtained on V (E) E = liom{'D*,E) by setting 
<j{x) = e o X o fj,, J{x) = Jc ° X and Vl — a o J . in particular, 

J{ri, ■■■ ,fn) = {Jen, ■■■ , JcTn), 

and if E is 3-dimensional, J{f\,- ■ ■ , r„) = {-^ Af\, - ■ ■ jt^tA r„). This defines 

1^1 1^1 

a (non-degenerate) hermitian structure on the subspace formed by the elements 
X = (rl,--- ,r^) e llom{V*,E) such that each belongs to the fixed space. 
In the sequel, whatever be a; e T) ® E, the set of all elements of the form 
\\x + A2i7(a;), Ai, A2 € M, will be called the complex line generated by x. 

As made clear by the very definition of C, the following hermitian structure 

plays an important role: 

Definition 3 The natural hermitian structure w^) ofV^ is defined by 

K^{u,v) = {iJ,{u},iJ,{v)), j{u,v) = {-v,u), uj^{u,v) = {-n{v),ii{u)). 

The Poisson structure The space EIom_|_ ((I?*)-^, I?^) of relative states is en- 
dowed with a Poisson structure whose symplectic leaves are the intersections of 
the submanifolds obtained by fixing the rank of £ and of those obtained by fix- 
ing the rotation invariants of the angular momentum. One fixes these invariants 
by fixing the traces of the iterates (of even order, those of odd order are equal 
to zero) of o f , which are equal to those of the iterates of C o e. 

3 Sundman's inequality 

3.1 The "components" of the angulcir momentum 

Recall that in endowed with its canonical basis, the vector product r^VAr 
by the vector V = (a, 6, c) is represented by the antisymmetric matrix 

/ -c b\ 
V = I c -a . 
\-b a / 

The scalar product V -W = {a,b,c) ■ {p, q, r) = ap + bq + cr can then be written 

V.W= ]^trace{V*W). 

Recalling that the natural coupling beween Hom{E,F) and Hom{E* ,F*) is 
given hy < f,ip >= trace{f o *<^) this motivates the following 

Definition 4 Let C : E* ^ E and Q : E ^ E* be antisymmetric (think 
of C as an angular momentum and Q as an instantaneous rotation). We call 
^{C,n) = ^trace{Co\l) the "component" of C "along" n. 
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Taking O = fic, we get 

Definition 5 The norm \C\ of the bivector C is defined by the formula 

\C\ = ^{C,nc) = \trace^-(CotY. 

When n = 2 or n = 3, |C| = ||C|| is the norm of the angular momentum vector 
(notice that Oc being normalized, corresponds to a vector of length 1). 
In the general case, if dimF = 2d, we have 

\C\ = \uii\ + \u)2\ H h l^dl, 

where the eigenvalues of (C o e)|ir are ±iwj, i = 1, • • • ,d. 

3.2 Complex Schwarz inequalities, Sundman's inequality 

Let (e, i7, = e o J) be a possibly degenerate Hermitian structure on the eu- 
clidean space {E,€). This means that for any f £ E we have the inequality 
||i7r||e < ||r||e instead of having the equality. 

As there is no risk of ambiguity, we denote by the same letter J the complex 
structure onV^E defined by the composition x ^ J ox (also denoted by Jx). 

Recalling the notations I = x - x, J = x - y, K = y ■ y, we have 

Proposition 6 For any z = (x, y) G 'D'^(E)E and any fl € Hom{E, E*) defining 
with € a degenerate Hermitian structure, one has 

IK-J^-^{C,flf>0. 

Equality occurs if and only if x belongs to the image of J and y belongs to the 

complex line generated by x. 

Proof. One computes 

(C, fi) = trace{—y o ^o^x o eo J) + trace{x o jjLo^y o eo J) = 2Jx ■ y. 

One concludes by the complex Schwarz inequaliity 

\\xW'\\yW'-{x-yf-{Jx-yf>0, 

which one proves by projecting y on the complex line generated by x and using 
the inequality WJxW < \\x\\, which is an equality only if x S ImJ) . It follows 
that equality occurs if and only if = ||a;|| and y belongs to the complex 

line generated by x. Choosing O = fic gives the 

Proposition 7 (Sundman's inequality) For any z = {x,y) G V^, one has 

IK-J^- |C|^ > 0. 

Equality occurs if and only if x belongs to the fixed space F and y belongs to the 
complex line generated by x. 
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3.3 Sundman's inequality as the best complex Schwcirz 
inequality 

Proposition 8 //, for some z = {x,y) £ V"^ E and some € Hom{E,E*) 
defining with e a degenerate Hermitian structure, equality occurs in the complex 
Schwarz inequality, the motion with initial condition z takes place within the 
fixed space F, on which the possibly degenerate hermitian structure J = oQ, 
coincides with Jc- 

Proof. The hypothesis imphes that there exist A, S K such that y = \x+vjx 
and that the /x e norm of Jx is the same as that of x\ it follows that ^7 is a 

non-degenerate hermitian structure on the image of z. 

From now on, we suppose that E — Imz is the space of motion. 

Without changing C = v{—ho*J + Joh) (where 6 = .x o /i o *a; is the inertia 

form), we may replace z = {x,y) by zo = (.xo,yo) {^/I'x, y/vj'x). The virtue 

of zq is to be a complex mapping from {T>*Y to E, where the complex structure 

on {V*fisj*='j-^:{^-r,)^{-r,,0-- 

Jozo = ZQOj*. 

Hence 

C = Zq O UJ^ O % = Zo O Kf, O j O % = -Zq O K^O % O *J = -{%)* K^J . 

As, by our assumption on E, *zq is injectivc, this implies that 

is symmetric and positive definite. In particular, C o J7, hence C, is an isomor- 
phism and the space of motion E coincides with the fixed space F. Moreover 

(C o *J o ef = {Co*Joe)o {-J o C o e) = -(C o ef . 

and hence 

-CoQ = Co*Q = Co*Joe = v/-(Coe)2 = -C o VLq, 
which implies that O = fic- 

3.4 Sundman's inequality and the decomposition of veloc- 
ities 

We give now an interpretation of Sundman's inequality which sheds light on the 
cases of equality: a configuration x E V (E) E being fixed, the space of velocities 
decomposes into the /x ® e-orthogonal sum of three subspaces, respectively 
(i) homothetic velocities (proportional to x), corresponding to changes of size, 
(a) rotational velocities, corresponding to solid body motions, 
(Hi) deformation velocities, corresponding to deformations of the shape: 

y = yh + yr + yd&'C"SiE. 
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This is the so-called Saari decomposition of velocity. 

One checks that = J^/I. Morever, as we have no a priori knowledge of 

the changes which affect the shape of the configuration we can only assert that 
llydiP ^ 0. Hence Sundnian's inequality is equivalent to ||?yr|P > |Cp//. and 
equality occurs if on the one hand there is no deformation velocity: yd = 0, and 
on the other hand = \C\^/I, which is equivalent to x and y belonging to a 
given complex line in the fixed space. 

Hence if a motion t i-^ x{t) = (fi(f),--- ,r„(i)) is such that the equality is 
verified at each instant t, there exists a fonction t X{t) G C, necessarily of 
class C°° (project one body on some axis), such that for all t, x{t) = X{t)x(t), 
where the multiplication is given by the hermitian structure defined by the 
angular momentum. The unicity of the solutions of differential equations implies 
that there exists a function t n- C,{t) G C, of class C°°, such that x{t) = C,{t)xQ. 
Using the invariance of U under rotation and its homogeneity, we get x = (^xq = 
Vf/(Ca;o) = j^VU{xo). Moreover, taking the scalar product with xq and using 



Euler identity xq ■ VU{xo) = -U{xo), we get 4|CP = 



CI 



_ Ujxo) 



Finally, we have 



VU{xo) = - — — xo, 
I{xo) 

^ A^o)|CI^' 



Definition 9 A configuration x is central if, when released with initial veloc- 
ities, it collapses homothetically on its center of mass. Equivalently, there exists 
A e M such that VU{x) = Xx (necessarily X = —j^), that is x is a critical 
point of U among configurations with fixed I. 

The first equation asserts that xq is a central configuration. In the case of three 
bodies, these configurations were described by Euler and Lagrange (see below) 
but describing the possible central configurations with more than three bodies 
is a very difficult problem. 

The second one is Kepler equation in C = M^. It implies that the complex 
homothetic motions x{t) = (^{t)xQ are Kepler-like: each body follows a Keplerian 
trajectory and, in particular, obeys Kepler laws. Notice that in the case of a 
(real) homothetic motion as in the above definition of central configurations, 
Sundman's equality becomes IK — = 0. 

Finally, it is not astonishing that the "simplest" (in the sense that Sundman's 
equality is realized) solutions of the n-body problem are exactly those which 
generalize the solutions of the 2-body problem or, what is the same, the solutions 
of the Kepler problem of attraction by a fixed center with the inverse square law, 
where there is no "shape" and the configuration stays constant up to similarity. 
In the next section, we briefly recall the integration of the Kepler problem in 
the case of negative energy. 
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3.5 Elliptic Kepler motions 

The Kepler equation in C = is 

■|ci 

We set 

C = re'\ I = Id' = r^ J = ReCC =U = rf,C = ImZt, K = |C|^ 

U=^=kp,H=^K-U = --, + 2Hc^ = fc2 - — = ee^. 
r 2 2a a 

As expected we have IK — — = 0, which we write in two ways 

aJ'^ + ir-kaf-k^a^e^ =0 or f"^ = 2H + 2kp - . 

The first equations defines an ellipse in the (r, J = rf) plane. On the other 
hand, if C ^ 0, the second equation defines an ellipse in the {r,p = ^) plane. 
Playing with the parametrization of these ellipses, we get the well-known (?) 
equations of elliptic (i.e. with negative energy) Kepler motion: 

c2 , 

r = —7- r, £ = ka(cosu — e), ri = kayl — e^sinu, 

A;(l + ecosu) 

t — to = kai {u — esinu) = ka^l, 

where ( = re'" = ^ + ir], and u, v, I are respectively the eccentric anomaly, the 
true anomaly and the mean anomaly. The motion takes place on an ellipse 
of eccentricity e and semi-major axis ka, the extreme cases e = and e = 1 
corresponding respectively to circular and coUinear motions. 



4 Simple motions originating from the symme- 
tries of the n-body problem 

An n-body solution for which Sundman's inequality is at each moment an equal- 
ity is such that the configuration of the bodies stays constant up to similarity. 
We shall prove a partial converse of this property: an n-body solution whose 
configuration stays constant up to similarity realizes equality in Sundman's in- 
equality except maybe in the case where it is a rigid motion in a space of 
dimension at least 4, in which case IK — ~ I^P is constant but maybe strictly 
positive. 



4.1 Rigid motions: balanced configurations 

A solution x{t) of the n-body problem is called rigid if the relative configuration 
Pit) = *x{t) o eo x{t) does not depend on time. The following proposition is a 
pleasant application of the form that we gave to the reduced equations (in fact 
I do not know of any other proof) : 
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Proposition 10 A motion is rigid if and only if it is a motion of relative 
equilibrium, that is if it defines an equilibrium of the relative equations. In 
other words, a solution of the n-body problem satisfies 7 = i/ and only if it 
satisfies £ = 0, that is ^ = 0, S + *A o P + P o A = 0, [A, p) = [A, P) = 0. 

Proof. One writes $ = P = P = P = and one notices that [A, [A, P)) = 
implies [A, P) = 0. 

This last equation, where only on the shape of the configuration and the forces 
comes into play, is important enough to deserve a name: 

Definition 11 An n-body configuration such that [A,P) = is called "balanced" 
(in french "equilibree"). 

The name comes from the fact that these configurations are exactly the ones 
which admit a relative equilibrium motion in a space of large enough dimension 
(2n — 2 is of course sufficient for n bodies): the attracting forces can be exactly 
balanced by the centrifugal forces. Indeed, given P > such that [A, p) = 0, one 
obtains a solution £ of the equation £ = hy setting j = p = 0, S = —2P o A. 
Notice that the positivity of S, which is equivalent to the negativity of the 
bilinear form A o fi = dU{P) on the image of P, is satisfied because of the 
attractivity of the Newton force. It follows that the quadratic form £ defined 
in this way is positive and hence of the form *z o e o z foTC z = {x,y) : V* E 
with dim E = 2 rank /?. 

The following proposition summarizes the properties of relative equilibrium 
motions. 

Proposition 12 A relative equilibrium motion always takes place in a space of 

even dimension and its configuration is balanced. It is a uniform quasi-periodic 
rotation of the absolute state z, that is: there exists a bivector fl € f\^E* , 
independant of time, such that z = oVto z. 

Proof. The equation 5 = of relative equilibria amounts to the antisymmetry 
of the tensor hoeo z, that is 

*2;oeoi e A^P^. 

As the motion takes place in the Image oi z : V'^ ^ E, we may suppose that 
z is surjective. Now the bivector ho to z has its support in Imh and we can 
define f2 e a'^E* as its image under that is 

^oVLo z = *zo eo z, which implies i = oQ,o z. 

The following diagramm summarizes the situation: 

{V*f > {V*f/Kerz — ^ E 

oeoz *jso€oi n 

< /m% ( *^ E* 
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We have yet to check 

1) that n does not depend on time, but this is a consequence of the unicity of 
solutions of differential equations (one can also directly check that f2 = from 
the fact that *z o o z = *z o e o i, being expressed with A, /3. 7, (5, p, is constant.) 

2) that the dimension of E is even, but this is a consequence of the non- 
degeneracy of fl on E (recall that we supposed that z is surjective). Indeed, if 
r E kerfl, then (*i o e)r = — (*z o n)r = 0. As i = {y. 2xA), this implies that 
er belongs to the intersection of the kernel of *y and the kernel of ^ o which 
is the same as the kernel of *a; because *A is non degenerate on the image of *x. 
Finally, ef € ker*z = 0. 

Other characterizations of balanced configurations. 

1) A relative configuration Pq is balanced if and only if it is a critical point of 
U among the configurations /3 such that h = fjbo (3 belongs to one and the same 
orbit of the democracy group 0{'D*), that is among the configurations /3 such 
that B = jjLo j3 has a given spectrum. Compare to the similar characterization 
of central configurations where only the trace I of B = fj, o j3 is fixed. 

2) A configuration x is balanced if and only if there exists a symmetric ho- 
momorphism T, : E ^ E* such that S/U{x) = o E o x. Indeed, recalling that 
V[/(x) = 2xo A and i^o*A = Ao jj, 

(i) if2xoA = oT,ox, we have 2BoA = fj.o^xo'Eox = 2fio*Ao*xoeox = 
2A o/io*a;oeoa: = 2A o B; 

(a) conversely, \i Ao B = B o Awe have found O such that z = o Q, o z, 
which implies WU{x) = (e~^ o O)^ o x and we set o S = (e~^ o 

Notice that, using this definition, the existence of relative equilibrium motion in 
a space whose dimension is twice the rank of x is even more transparent: after 
noticing that one can chose E such that o S preserves the image of x, it 
remains only to complexify this image in order to find a space in which e~"^ o S 
has an antisymmetric square root o f2. 

4.2 Homographic motions: central configurations 

One calls homographic a solution z{t) = [x{t),y{t)) of equations (N) such that 
there exists a real function of the time and a relative configuration /?o with 
P{t) = z^(i)^/3o, where /3{t) — *x{t) o e o x{t). They comprise two important 
particular cases : 

- (real-)honiothetic solutions solutions z{t) = [x{t),y{t)) such that there ex- 
ists a real function 1/ of time and an absolute configuration xq with x{t) = i^{t)xo- 

-rigid solutions already studied, such that the relative configuration (3{t) does 
not depend on time. 

We shall see that, at least when they are not relative equilibria, they are exactly 
the Kepler-like motions which we defined above as the ones which realize the 

equality in Sundman's inequality and hence are complex homothetic. More pre- 
cisely, the setting that we introduced allows giving short proofs of what follows. 
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Proposition 13 The configuration of a non rigid homographic motion is cen- 
tral. 

This proposition becomes false if the Newton potential is replaced by the Jacobi- 
Banachiewitz potential. 

The definition of a central configuration was given when studying the cases of 
equality in Sundman's inequality. Notice that a central configuration is balanced 
because one deduces from the definition that 2f3 o A — X/3 ^ *(A/3) = 2*^ o /?. 

The description of (real-)homothetic solutions is easy : for any normalized (i.e. 
7 = 1) central configuration xq and any real solution of the differential 
equation ( = 2kU{xq)\(\'^'^~'^C (essentially (N) in the one-dimensional case), 
x(t) = Ci't)x is a homothetic solution. Moreover, every homothetic solution is 
of this type. 

The following proposition shows that non homothetic motions whose configura- 
tion is central (this is automatic except for relative equilibria in dimension at 
least 4) are indeed complex homothetic : 

Proposition 14 The space of motion Im z of a homographic, non homothetic, 
solution with central configuration, coincides with the fixed space. For the com- 
plex structure on Hom(I?*,Imz) induced by the angular momentum, y is at any 
time a complex multiple of x : if Xq — ||a;(0)||^"'^a;(0) is the normalized initial 
configuration, x(t) — Cit)xo where C, is a complex function of the time satisfying 
^ = — t/(xo)|C|~^C- Inversely, any complex solution of this differential equation 
gives rise to a complex homothetic solution. 

One deduces from this proposition that a homographic motion with central 
configuration, in particular any non rigid homographic motion, is a generalized 
Keplerian motion: all bodies describe similar conies around the center of mass. 

Proof. [Proof of proposition [13] Let z{t) be a solution whose relative configu- 
ration I3{t) is of the form /?(t) — i/{t)'^Po. We suppose that the initial relative 
configuration jSg is normalized in such a way that Iq = tracer o /3 = 1. Then 
I{t) = vitY and by homogeneity, we have at any time 

We choose a /x^-'^-orthonormal basis {ui, ■ ■ ■ of the vector space T)* in 

which the ^-symmetric endomorphism Aq is represented by a diagonal matrix 
Aq = diag{ai, • • • , a„_i). We call foi, ■ ■ • , &n-i the diagonal terms of the matrix 
representing Bq = ^ o f3o in this basis. Recalling that /3 = x*e, we have 

lk(w»)||? = \\u^\\l = I\\u,\\l^ = Ih, 

the last equality coming from the fact that the basis {ui, • • • , is orthonor- 

mal. Now, for each index i such that bi ^ 0, we define 

_ 1 

Si = 6, ^x{ui) e E. 
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The functions t H> Si{t) satisfy the equations 

si = 2\\s^\\~^aiSi and ||si||^ = /. 

In other words, they are solutions of a priori different Kepler equations but their 
norms coincide at each instant. If the Si(t) are circular solutions, nothing else 
can be concluded. But if they are not circular, that is if I is not constant, this 
implies that the Kepler equations are all the same: = a is independant of i. 
Indeed, if we set 

Ii = Ji =< Si, si >, Ki = ||sl|p. Hi = ^Ki + 2ail~^, 

the "energies" Hi are constant. Now, by hypothesis, all the /, coincide with the 
function / and hence all the Ji coincide as well as all the Ji. By the Lagrange- 

Jacobi identity, we have Ji — 2Hi — iail,^ ^ , hence 

2{H, - Hj) = j,- jj + 2a,I~^ - = 2{a, - a^)!'^ 

can be constant only if either / is constant or = aj. This proves that, if I is 
not constant, the restriction of A to Im x is of the form Alimx = aid and hence 
that a; is a central configuration. Note that the above computation amounts to 
saying the obvious fact that if an elliptic Keplerian motion is not circular, the 
function I determines its period. Here also, the conclusion does not hold for 
the Jacobi-Babachiewitz case of a force proportional to the inverse cube of the 
distance. 



Proof. [Proof of proposition 14 We consider now a homographic motion 
which is not homothetic and whose configuration is central (recall that this last 
condition is automatically satisfied if the motion is not a relative equilibrium) . 
Setting I — J =< x,x >, K = the hypothesis implies IK — J"^ ^ 0. 

We have 

m = I{t)Po, lit) = J{t)(3o, Sit) = Kit)l3o. 
To prove the last identity, we write 

27 = //3o = 26 + 21-^% o /3o + /?o o Ao). 
By hypothesis, 2a;o°^o = 2\xq, hence /SqoAo = A/3o = A*/3o = *Ao o /Sg. Finally, 

\Jpo + p Kfio 

We want to prove that the velocity configuration y is a each time the composition 
y = (^x of the configuration x with the multiplication by some complex number 
( = ^ + iri in some hermitian structure onV^E. As ix must be /i (g) e-orthogonal 
to X, one must have 

i= J, irix = y~ jx. 
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Using the fact that the multipHcation by i in 2? (g) is antisymmetric, one com- 
putes hjey = {ij- + Ir]^)j3Q, which impHes that (choosing the complex structure 
in such a way that 77 > 0) , 



y 



1 r 

7 



Setting c — \J IK — ^ 0, we normahze z by replacing it by zq — {xo,yo) in 
such a way that yo = ixo by setting 



Xo 



Lemma 15 There exists a hermitian structure Jo on E such that 

yo = JQ° xq. 

Proof. Let zq = (a;o,yo)- We are looking for a hermitian structure on E, 
i.e. an e-isometry Jo whose square is —Ids, which defines the multiplication 
by j = V— 1. It is necessarily of the form Jo = e^^ o rip, where Vlo : E ^ E* 
is antisymmetric. Let j* : (V*)^ (2?*)^ be the natural complex stucture 
defined in section 2.4 by jo{£,,v) — {~V^0- If ^ Jo ° xo, we must have 
zo° i* ~ {yo, ~xo) = Jo° zo, hence £0° j* = *zo o 51o o zq must be antisymmetric. 
This is indeed the case: 

-P -c/3| 



c/3o -P 



£o°f = zooeozoof 



Hence, as in proposition |12| supposing zg surjective, one deduces the existence 

oi ^lo '■ E E* antisymmetric such that zq o j* = e^^ o JIq o zq. 

Finally, yo — e^^ o JIq ° a;o and — zq = zq o (j*)^ — (e^^ o flo)'^ o zq which implies 

that (e~^ o r^o)^ = —Ids because zq is surjective. 

We summarize the situation in the following commutative diagram: 



(2?*)2 ^ 


E 


^ E 


1'- 








E 





Coming back to the original configuration z = {x,y), we have proved that 
y = J [J + ic] X, where the action of C on the configuration space is given by 
the complex structure Jo = e^^ o ilo. 

The angular momentum of z = {x, y) = (x, j{Jx + ce^^ o fig o x)) is 
C = j{b o ilo o e^^ + e^^ o ilg o 6), where 6 = x o /i o *x (= inertia = S o e^^). 
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As the inertia endomorphism S ^ b o e has the same trace / as /i o /3, we see 
that the complex Schwarz inequality for f^o is an equality : 

11 c / 

y ■ {Jo ° 2;) - (C, f^o) — -trace('C o Qq) — —traceS = c = V IK — J^. 

From proposition |8] we deduce that fio = i^c • The end of the proof of Proposi- 
tion [14] follows from section |331 



4.3 Equations for the balanced and central configurations 

To write down usable equations of these configurations it is convenient to rep- 
resent bilinear forms on I? or on by n x n matrices. Such a representation is 
unique in the case of T) but ambiguous in the case of V* where we have to chose 
an extension of the form to R". Representing /3 by the matrix of general term 
— c?J7 by the matrix of general term — and /i^^ by the diago- 

nal matrix of the m~^, one finds that 11 = /3 o dU o pT^ = /3 o ^4 e Hom(2?*, I?) 
is represented by the matrix whose coefficients are the 

To get convenient coordinates for the antisymmetric part of 11, one notices that 
the exterior product by (1, • • • , 1) factorizes through an embedding of /\^ P in 
A M". The coordinates of the image of the bivector 11 — *n by this injection 
are the 

Pi/ik Pii ^ ^ Pki Pih Phi Pi', 



9 ^ijfei 



and 



with i < j < k, where 

1 /rrii 1 1 rrij 1 / rrik 

^ ijk — det I Sjk Ski Sij Ski Sij Sjk Sij Sjk Ski 

dU/dsjk dU /dski dU /ds,j 

1 1 1 

y.' = det I Sjk + Sii Skt + Sji s^j + Ski 
' du_ J_9(L 1 au 

dsn rrij dsji nik dski 

The equations Pijk = 0, i < j < k, define the balanced configurations. Of 
course, if n is strictly greater than three 3, we get too many equations but this 
is nevertheless the best way of writing down the equations. 

As U — ^i<:i^j<:„'mi'mj^{sij), where $ is defined by $(s) = Qs~^^^ in the 
Newtonian case, the above equations are linear in the masses ! A most important 
property of $ is the concavity of its derivative if. For three bodies, it implies 
immediately the following proposition : 
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Proposition 16 In the Newtonian case, a configuration of three bodies of equal 
masses is balanced if and only if it is isosceles. 

Figure 1 shows the shape sphere, that is the space (IR^)^/50(2) of triangles 
of constant moment of inertia, up to oriented similarity in (see |C1| 1. The 
level lines of the potential in the case equal masses are indicated; in this rep- 
resentation, the level lines of the area function are horizontal circles. Using 
the symmetries of the potential function and the variational characterization of 
balanced configurations given above (critical points of U among triangles with 



moment of inertia and area fixed), this gives another proof of Proposition 16 




Figure 1: the shape sphere 



5 Simple motions in odd dimensions related to 
relative equilibria 

All the motions studied up to now take place in a space of even dimension. 
In odd dimensions it turns out that a good measure of the "simplicity" of a 
periodic solution is its Lagrangian action. Moreover, in some cases, periodic 
solutions which minimize the action among loops of configurations sharing a 
certain type of symmetry are naturally related to relative equilibria through 
families of quasi-periodic solutions. 

5.1 The Hip-Hop in 

What are the "simplest" periodic solutions of the 4-body problem ? In or 
M.^, it seems natural to look for a candidate among the relative equilibria of a 
central configuration, as they turn Sundman's inequality into an equality. In 
order to distinguish between the various central configurations, one can compute 
the Lagrangian action Jq (y + U)dt for a given period T and notice that it is 
minimum when the configuration is a minimum of the potential U restricted to 
the configurations with a given moment of inertia /. 
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In this selects the square and in K'' the regular tetrahedron. But what 
about motions in K'^ ? In |CV| a solution called the Hip-Hop is proved to 
exist which, in some sense, interpolates between these two. It minimizes the 
Lagrangian action among loops of configurations with a given period T which 
are invariant under some action of the group Z2 x Z4 (or Z2 x D4) . The projection 
of the configuration on some plane (say the horizontal plane) is always a square, 
whose size varies gently, while the two diagonals make one vertical oscillation per 
period in opposite directions. In particular, twice per period the configuration 
is a square and twice per period it is a regular tetrahedron. 

One can reasonably conjecture that imposing the sole italian symmetry, that 
is the identity x{t — ^) = —x{t) for all t, would give rise to the same solutions. 
In any case, it can be proved that such minimizers have no collisions (i.e. that 
they are true solutions) and that, for the spatial problem, they are never planar 
; moreover, this is true for an arbitrary number > 4 of bodies and any system 
TOi, • • • , mjv of positive masses [C4j . 



5.2 From to M'^: the unchained polygons 

The equal mass Hip-Hop is naturally related to the relative equilibrium of the 
square through a family of quasi-periodic (in fact periodic modulo rotation, 
i.e. periodic in rotating frames) solutions in M.^ (see |TV[ICF| ): in a well chosen 
family of rotating frames, the members of the family look like Hip-Hop solutions 
parametrized by the amplitude of their vertical oscillations. 

More generally, the (planar) relative equilbria of N equal masses forming a 
regular iV-gon are related in the same way to planar periodic solutions which 
can be either generalized Hip-Hops or choreographies |CF| . 




Figure 2: the original Hip-Hop {p = 2) 
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Figure 3: the Hip-Hop family (all frames are rotating around the vertical axis 

except the one at right) 

5.3 Hip-Hops in higher dimensions 

In higher dimensions, the situation seems to be exactly similar : given 2p bodies, 
the action minimizers with Italian symmetry constraint are the relative equilibria 
of the regular 2p-simplex in R^'' and one can conjecture that: 

(i) h\ M^P"^ they are relative equilibria of a central configuration made of 
a copy of the regular p-simplex in each one of two orthogonal subspaces of 
dimension p — 1; 

(ii) in M^P^^ they are of Hip-Hop type, with the two diagonal segments of 
the p — 2 case replaced by two regular p-simplices which would have opposite 
periodic motions along the K factor of K^^"^ = R'^p-^ x M. 

5.4 From one Hip-Hop to the other through the fourth 
dimension 

Other equal mass Hip-Hop solutions exist which have different symmetries (see 
|CVI IC4I ICF| ). Figure 4 shows how two of them are simply related through 
families of quasi periodic solutions in M^. The key fact is that, while the reg- 
ular tetrahedron is a central configuration, Z/4Z-symmetric tetrahedra as well 
as Z/3Z-symmetric tetrahedra are balanced configurations which admit quasi- 
periodic relative equilibrium motions in (see [C5] , this means some degen- 
eracy occurs, as a general balanced configuration of 4 bodies would allow only 
quasi-periodic relative equilibrium solutions in M^). 

1) the Hip-Hop family of quasi-periodic solutions in M.^ described in figure 3 
leads from the Hip-Hop to the planar periodic relative equilibrium of the square; 

2) a family of quasi-periodic relative equiliibria in of Z/4Z-symmetric less 
and less fiat tetrahedra links the relative equilibrium solution of the square to 
a periodic relative equilibrium solution of the regular tetrahedron; 

3) a family of periodic relative equilibria in of the regular tetrahedron, 
corresponding to a family of complex structures in M**, leads to a new periodic 
relative equilibrium solution of the regular tetrahedron; 

4) a family of quasi-periodic relative equiliibria in of Z/3Z-symmetric 
more and more fiat tetrahedra links the new relative equilibrium solution of the 
regular tetrahedron to the periodic planar relative equilibrium solution of the 
equilateral triangle with the fourth mass at its center; 

5) a family of quasi-periodic solutions in of Hip-Hop type with Z/3Z- 
symmetry ends with a periodic Hip-Hop solution with Z/3Z-symmetry. 
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The eigenvalues that the angular momentum of a relative equilibrium can take 
and the ones at which a periodic relative equilibrium with central configuration 
may undergo a bifurcation to a family of quasi-periodic relative equilibria with 
balanced configurations are studied in |C2| . The first question amounts to a 
purely algebraic problem related to Horn's problem, which is now completely 
understood |CJj . 




Figure 4: from a Z/4Z-symmetric Hip-Hop to a Z/3Z-symmetric one 
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